11 Lorenz equations 



In this lecture we derive the Lorenz equations, and study their behavior. 

The equations were first derived by writing a severe, low-order truncation of 
the equations of R-B convection. 

One motivation was to demonstrate SIC for weather systems, and thus point 
out the impossibility of accurate long-range predictions. 

Our derivation emphasizes a simple physical setting to which the Lorenz 
equations apply, rather than the mathematics of the low-order truncation. 

See Strogatz, Ch. 9, for a slightly different view. This lecture derives from Tritton, Physical Fluid 
Dynamics, 2nd ed. The derivation is originally due to Malkus and Howard. 

11.1 Physical problem and parameterization 

We consider convection in a vertical loop or torus, i.e., an empty circular 
tube: 

cold 




hot 



We expect the following possible flows: 

• Stable pure conduction (no fluid motion) 

• Steady circulation 

• Instabilities (unsteady circulation) 
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The precise setup of the loop: 
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Tq+Tj (external) 
4> = position round the loop. 

External temperature Te varies linearly with height: 

T E = T - T x z/a = T + T 1 cos , 
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(24) 



Let a be the radius of the loop. Assume that the tube's inner radius is much 
smaller than a. 

Quantities inside the tube are averaged cross-sectionally: 

velocity = q = q((f), t) 
temperature = T = T(0, t) (inside the loop) 



As in the Rayleigh-Benard problem, we employ the Boussinesq approximation 
(here, roughly like incompressiblity) and therefore assume 

Thus mass conservation, which would give V • u in the full problem, here 
gives 
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Thus motions inside the loop are equivalent to a kind of solid-body rotation, 
such that 

q = q(t). 

The temperature T(0) could in reality vary with much complexity. Here we 
assume it depends on only two parameters, T 2 and T3, such that 

T - T = T 2 cos (J) + T 3 sin 0. (26) 

Thus the temperature difference is 

• 2T2 between the top and bottom, and 

• 2T3 between sides at mid-height. 

T2 and T3 vary with time: 

T 2 = T 2 (t), T 3 = T 3 (t) 

11.2 Equations of motion 
11.2.1 Momentum equation 

Recall the Navier-Stokes equation for convection: 

— * Y —* 

— — h u ■ Vu = — Vp — gaAT + uV 2 u 
at p 

We write the equivalent equation for the loop as 

dq 1 dp 

at pa oq> 

The terms have the following interpretation: 

• u — > q 

• u • Vm — > since dq/d(j) = 0. 

Ill 



(27) 



• Vp — > by transformation to polar coordinates. 

• A factor of sin0 modifies the buoyancy force F = ga(T 
the tangential component: 



Tq) to obtain 




Fsin§ 



The sign is chosen so that hot fluid rises. 

• T is a generalized friction coefficient, corresponding to viscous resistance 
proportional to velocity. 

Now substitute the expression for T — Tq (equation (26)) into the momentum 
equation (27): 



dq 1 dp 

dt pa d(j) 



+ ga(T2 cos (f) + T3 sin 0) sin <j) — Tq 



Integrate once round the loop, with respect to 0, to eliminate the pressure 
term: 

da f 2lT 

2tt-± = ga / (T 2 cos sin + T 3 sin 2 0)d0 - 2nYq. 
dt Jo 

The pressure term vanished because 

Jo d0 

i.e., there is no net pressure gradient around the loop. 



The integrals are easily evaluated: 

r 2n 



1 



/ cos (j) sin <j) d(j) = — sin <j) 
'0 2 



2tt 



= 



and 



sin 0d0 = 7r. 
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Then, after dividing by 27r, the momentum equation is 

|-r 9 + Sf (28) 

where we have written dq/dt instead of dq/dt since dq/d(f) = 0. 

We see that the motion is driven by the horizontal temperature difference, 
2T 3 . 



11.2.2 Temperature equation 

We now seek an equation for changes in the temperature T. The full tem- 
perature equation for convection is 

dT 

~ +u-VT = kV 2 T 



at 

where k is the heat diffusivity. 

We approximate the temperature equation by considering only cross-sectional 
averages within the loop: 

Here we have made the following assumptions: 

• RHS assumes that heat is transferred through the walls at rate 

K (-^external ^internal)- 

• Conduction round the loop is negligible (i.e., no V 2 T). 

• 1 1? is the product of averages, not (as it should be) the average of a 
product; i.e., q is taken to be uncorrected to dT/dcf). 

Recall that we parameterized the internal temperature with two time-dependent 
variables, T 2 (t) and T^(t). We also have the external temperature Te varying 
linearly with height. Specifically: 

T E = r o + r lC os0 
T — To = T2 cos + T3 sin 
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Subtracting the second from the first, 

T E -T =(Ti-T 2 ) cos0-T 3 sin0. 
Substitute this into the temperature equation (29): 

dT dT 

— — - cos (f> H — — - sin (f> — -T 2 sin + -T 3 cos = tf(Ti - T 2 ) cos - AT 3 sin 0. 
dt dt a a 

Here the partial derivatives of T have become total derivatives since T 2 and 
T 3 vary only with time. 

Since the temperature equation must hold for all 0, we may separate sin^> 
terms and cos terms to obtain 

. dT 3 qT 2 

sm0: — = -KT 3 

dt a 

cos0: ^ + = K{ Tl -T 2 ) 
dt a 

These two equations, together with the momentum equation (28), are the 
three o.d.e.'s that govern the dynamics. 

We proceed to simplify by defining 

T 4 (t) = T 1 -T 2 (t), 

which is the difference between internal and external temperatures at the top 
and bottom — loosely speaking, the extent to which the system departs from 
a "conductive equilibrium." Substitution yields 



dt a a 

dT 4 T ^ qT 3 

-f = -ifT 4 + — 
dt a 



11.3 Dimensionless equations 



Define the nondimensional variables 



x = Q y = 9aTs gaT 4 

' ah" 1 2(iVI\" " ~~ 2aTK 
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Here 

X = dimensionless velocity 

Y = dimensionless temperature difference between up and down currents 

Z = dimensionless departure from conductive equilibrium 



Finally, define the dimensionless time 

t' = tK. 



Drop the prime on t to obtain 

dX 



dt 

dY 
~dt 

dZ 



= -PX + PY 
= -Y + rX-XZ 
= -Z + XY 



dt 

where the dimensionless parameters r and P are 

qolT\ . . . . ,, 

r = — — — = Kayleigh number 

r 

P = — = "Prandtl number" 
K 

These three equations are essentially the same as Lorenz's celebrated system, 
but with one difference. Lorenz's system contained a factor b in the last 
equation: 

^ = - bZ + XY 
dt 

The parameter b is related to the horizontal wavenumber of the convective 
motions. 



11.4 Stability 



We proceed to find the fixed points and evaluate their stability. For now, we 
remain with the loop equations (6 = 1). 
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The fixed points, or steady solutions, occur where 

X = Y = Z = 0. 

An obvious fixed point is 

X* = y* = Z* = 0, 

which corresponds, respectively, to a fluid at rest, pure conduction, and a 
temperature distribution consistent with conductive equilibrium. 

Another steady solution is 



x* = y* = ±Vr^i 

Z* = r-l 

This solution corresponds to flow around the loop at constant speed; the ± 
signs arise because the circulation can be in either sense. That sgn(X) = 
sgn(y) implies that hot fluid rises and cold fluid falls. 

Note that the second (convective) solution exists only for r > 1. Thus we 
see that, effectively, r = Ra/Ra c , i.e., the convective instability occurs when 
Ra > Ra c . 

As usual, we determine the stability of the steady-state solutions by deter- 
mining the sign of the eigenvalues of the Jacobian. 



Let 




Then the Jacobian matrix is 



d<f>j 



0* = 



-P +P 
r-Z* -1 
Y* x* 






-X" 

-1 



The eigenvalues a are found by equating the following determinant to zero: 

= 



-O + p) p o 

r-Z* -0+1) -X* 

y* x* -O + i) 
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For the steady state without circulation (X* = Y* = Z* = 0), we have 

= 0. 



O + p) p o 

r -0+1) 

o -O+i) 



This yields 



or 



-0 + P)0 + l) 2 + rPO + 1) = 
+ 1) [a 2 + a(P + 1) - P(r - 1)] = 0. 



There are three roots: 



°x = -1 

-(P + l) v / ( p + 1 ) 2 + 4p ( r " 1 ) 



2 2 
As usual, 

Re{o"i, a"2, and 0-3} < =4> stable 
Re{o"i, d2, or 0-3} > =^> unstable 

Therefore X* = Y* = Z* = is 

stable for < r < 1 
unstable for r > 1 

We now calculate the stability of the second fixed point, X* = ±y/r — 1, 
y* = ± v / f T7 T, Z* = r - 1. 

The eigenvalues a are now the solution of 





■O + P) P 

1 -O + i) s 
s s -O + i) 



= 0, 



(Explicitly, 

-0 + p)0 + 1) 2 - ^s 2 - s 2 + p) + ^0 + 1) 

+ 1) [<r 2 + a(P + 1)] + aS 2 + 2PS 2 





0.) 
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The characteristic equation is cubic: 

d 3 + a 2 {P + 2) + a{P + r) + 2P(r - 
This equation is of the form 

a 3 + Aa 2 + Ba + C = 
where A, 5, and C are all real and positive. 

Such an equation has either 

• 3 real roots; or 

• 1 real root and 2 complex conjugate roots, e.g., 



1) = 




or 



3 real roots 




(30) 



1 real and 2 complex conjugates 



Rearranging equation (30), 

a (a 2 + B) = -Aa 2 - C < 0. 

positive real negative real 

Consequently any real a < 0, and we need only consider the complex roots 
(since only they may yield Re{cr} > 0). 

Let a i be the (negative) real root, and let 

cr 2j 3 = a ± i(3. 

Then 

(a — (Ji)(a — a — i/3)(a — a + i/3) = 

and 

A = -((7i + 2a) 
B = 2aai+a 2 + (3 2 
C = -ct^ + Z? 2 ) 
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A little trick: 

C-AB = 2a [Oi + a) 2 + /3 2 ] . 

positive real 

Since a is the real part of both complex roots, we have 

sgn(Re{cr 2 , 3 }) = sgn(a) = sgn(C - AB). 
Thus instability occurs for C — AB > 0, or 

2P(r - 1) - (P + 2)(P + r) > 0, . 

Rearranging, 

r(2P - P - 2) > 2P + P(P + 2) 
and we find that instability occurs for 

P(P + 4) 

This condition, which exists only for P > 2, gives the critical value of r for 
which steady circulation becomes unstable. 

Loosely speaking, this is analogous to a transition to turbulence. 
Summary: The rest state, X* = Y* = Z* = 0, is 

stable for < r < 1 
unstable for r > 1. 

The convective state (steady circulation), X* = Y* = ±y/r — 1, 
Z* = r - 1, is 

stable for 1 < r < r c 
unstable for r > r c . 

What happens for r > r c ? 

Before addressing that interesting question, we first look at contraction of 
volumes in phase space. 
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11.5 Dissipation 



We now study the "full" equations, with the parameter 6, such that 

Z = -bZ + XY, b>0. 
The rate of volume contraction is given by the Lie derivative 
1 dV 



V dt ^ d(j>. 
For the Lorenz equations, 



= E*X' * = 1,2,3, 1 = X,0 2 = y,0 3 = ^. 



dx dY dz „ , , 
ox + W + oz = - p - l - b - 

Thus 

dV 

- = - {P + l + b) V 
which may be solved to yield 

V(t) = V(0)e^ p+1+b ^. 
The system is clearly dissipative, since P > and b > 0. 

The most common choice of parameters is that chosen by Lorenz 
P = 10 

b = 8/3 (corresponding to the first wavenumber to go unstable). 
For these parameters, 



V{t) = V{0)e~ 

41 R 

Thus after 1 time unit, volumes are reduced by a factor of e~~ ~ 10~ . The 
system is therefore highly dissipative. 



11.6 Numerical solutions 

For the full Lorenz system, instability of the convective state occurs for 

P(P + 3 + b) 



r > r r = 



P-l-b 
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For P=10, b=8/3, one has 

r c = 24.74. 
In the following examples, r = 28. 

Time series of the phase-space variables are shown in 

Tritton, Fig 24.2, p. 397 



• X(t) represents variation of velocity round the loop. 

— Oscillations around each fixed point X+ and X*_ represent variation 
in speed but the same direction. 

— Change in sign represents change in direction. 

• Y{t) represents the temperature difference between up and downggoing 
currents. Intuitively, we expect some correlation between X{t) and Y{t). 

• Z{t) represents the departure from conductive equilibrium. Intuitively, 
we may expect that pronounced maxima of Z (i.e., overheating) would 
foreshadow a change in sign of X and Y, i.e., a destabilization of the 
sense of rotation. 

Projection in the Z-Y plane, showing oscillations about the unstable connec- 
tive fixed points, and flips after maxima of Z: 

BPV, Fig. VI.12 

A 3-D perspective, the famous "butterfly:" 
BPV, Fig. VI.14 

Note the system is symmetric, being invariant under the transformation X — > 
-X, Y -> -Y, Z -> Z. 

A slice (i.e., a Poincare section) through the plane Z = r — 1, which contains 
the convective fixed points: 

BPV, Fig. VI.15 
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• The trajectories lie on roughly straight lines, indicating the attractor 
dimension d ~ 2. 

• These are really closely packed sheets, with (as we shall see) a fractal 
dimension of 2.06. 

• d ~ 2 results from the strong dissipation. 

Since d ~ 2, we can construct, as did Lorenz, the first return map 

Zk+l = f(Zk), 

where Zk is the kth maximum of Z(t). The result is 
BPV, Fig. VI.16 

(These points intersect the plane XY — bZ = 0, which corresponds to Z = 0.) 

The first-return map shows that the dynamics can be approximated by a 1-D 
map. It also reveals the stability properties of the fixed point Z = r — 1: 

BPV, Fig. VL17 

Finally, sensitivity to initial conditions is documented by 
BPV, Fig. VI.18 

11.7 Conclusion 

The Lorenz model shows us that the apparent unpredictability of turbulent 
fluid dynamics is deterministic. Why? 

Lorenz's system is much simpler than the Navier-Stokes equations, but it is 
essentially contained within them. 

Because the simpler system exhibits deterministic chaos, surely the Navier- 
Stokes equations contain sufficient complexity to do so also. 
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Thus any doubt concerning the deterministic foundation of turbulence, such 
as assuming that turbulence represents a failure of deterministic equations, 
is now removed. 

A striking conclusion is that only a few (here, three) degrees of freedom are 
required to exhibit this complexity. Previous explanations of transitions to 
turbulence (e.g., Landau) had invoked a successive introduction of a large 
number of degrees of freedom. 
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